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Abstract 

Quantum computing was so far mainly concerned with discrete problems. Recently, 
E. Novak and the author studied quantum algorithms for high dimensional integra- 
tion and dealt with the question, which advantages quantum computing can bring 
over classical deterministic or randomized methods for this type of problem. 

In this paper we give a short introduction to the bàsic ideas of quantum computing 
and survey recent results on high dimensional integration. We discuss connections 
to the Monte Cario methology and compare the optimal error rates of quantum 
algorithms to those of classical deterministic and randomized algorithms. 



Introduction 

One of the most challenging qüestions of today, in the overlap of computer sci- 
ence, mathematics, and physics, is the exploration of potential capabilities of 
quantum computers. Milestones were the algorithm of Shor [26], who showed 
that quantum computers could factor large integers efficiently (which is widely 
believed to be infeasible on classical computers) and the quantum search al- 
gorithm of Grover [10], which provides a quadratic speedup over deterministic 
and randomized classical algorithms of searching a database. 



So far research was mainly concentrated on discrete problems like the above 
and many others one encounters in computer science. Much less is known 
about computational problems of analysis, including such typical field of ap- 
plication of Monte Cario methods as high dimensional integration. We seek to 
understand how well these problems can be solved in the quantum model of 
computation (that is, on a - hypothetical - quantum computer) and how the 
outcome compares to the efficiency of deterministic or Monte Cario algorithms 
on a classical (i. e. non-quantum) computer. 
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Abrams and Williams [2] suggested first ideas about quantum integration al- 
gorithms. A systematic study was begun by Novak [23], who considered inte- 
gration of functions from Hòlder spaces. This line of research was continued by 
the author [14], where quantum algorithms for the integration of L p -functions 
and, as a key prerequisite, for the computation of the mean of p-summable 
sequences were constructed. In [14] a rigorous model of quantum computation 
for numerical problems was developed, as well. The case of integration of func- 
tions from Sobolev spaces is considered in [15], and more on the computation 
of the mean is presented in [17]. These papers also established matching lower 
bounds. A short survey of first results can be found in [16]. Path integration 
is studied by Traub and Wozniakowski in [30] . 

Combining these results with previous ones of information-based complexity 
theory about the best possible ways of solving the respective problems de- 
terministically or by Monte Cario on classical computers, we are now in a 
position to fairly well answer the question where quantum computation can 
provi de a speedup in high dimensional integration and where not. There are 
cases where quantum algorithms yield an exponential speedup over determin- 
istic algorithms and a quadratic speedup over randomized ones (on classical 
computers). 

Moreover, there is a close connection of quantum algorithms with Monte Cario: 
While computations are carried out on superpositions of classical states of 
qubit systems, and thus in high parallelism, the result can only be accessed 
through a measurement, which destroys the superposition and outputs any one 
of the superposed states - with a certain probability. Thus, these algorithms 
are probabilistic, Monte Cario, while, on the other hand, completely different 
laws govern the computation. Nevertheless various Monte Cario techniques 
can be put into use to construct quantum algorithms (which then, combined 
with special quantum techniques, outperform their classical counterparts) . 

We start with a brief introduction to the ideas of quantum computing. Then we 
consider the question, what quantum computers could do in numerical (Monte 
Cario related) problems and survey recent results on summation and high 
dimensional integration. We discuss how to use Monte Cario methodology for 
the development of quantum algorithms. Complexity issues and comparisons 
of the potential of deterministic, randomized and quantum algorithms are 
considered, as well. 

For further reading on quantum computation we recommend the surveys by 
Aharonov [1], Ekert, Hayden, and Inamori [8], Shor [28], and the monographs 
by Pittenger [24], Gruska [12], and Nielsen and Chuang [21]. For notions and 
results in information-based complexity theory see the monographs by Traub, 
Wasilkowski, and Wozniakowski [31] and Novak [22], and the survey Heinrich 
[13] of the randomized setting. 
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1 A Short Introduction to Quantum Computing 



1.1 History 



The first ideas of using quantum devices for computation were expressed at 
the beginning of the eighties by Manin [18], see also [19], and Feynman [9]. 
They observed that simulating quantum mechanics on a classical computer 
is extremly hard, probably infeasible, since it leads to differential equations 
whose dimensions are exponential in the number of system components. To 
overcome this, they suggested the idea to simulate quantum mechanics using 
quantum devices itself. 

In 1985, Deutsch [7] developed the rigorous theoretical basis of quantum com- 
putation - the model of a quantum Turing machine, which became the so far 
most serious and still standing challenge to the Turing-Church Thesis (the 
latter stating that, roughly, every reasonable physical computing device can 
be simulated with only polynomial increase of resources on a classical Turing 
machine). 

A breakthrough for quantum computing happened in 1994, when Shor [26] 
showed that efficient factorization of integers would be possible on a quantum 
computer, which in turn, would mean the possibility of breaking the foremost 
públic key codes like the RSA cryptosystem. Another fundamental contribu- 
tion was Grover's [10] discovery of an efficient quantum search algorithm in 
1996 (we give some further comments on both Shor's and Grover's algorithms 
later on). 

Since then we witness an explosion of efforts, broad research on quantum algo- 
rithms for all kinds of (mostly discrete) problems, on quantum cryptography, 
and quantum information theory. Physicists are intensively working on how to 
construct quantum computers, that means, finding quantum mechanical Sys- 
tems that can be manipulated to fulfill the abstractly proposed requirements. 
Systems with a few qubits are already successfully realized in laboratories. 

Important forerunners for the development of quantum algorithms for Monte 
Cario related, numerical problems were the work of Boyer, Brassard, H0yer, 
Mosca, and Tapp [4], [5] on quantum counting and the results of Beals, 
Buhrman, Cleve, Mosca and de Wolf [3] and Nayak and Wu [20] on lower 
bounds. 
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1.2 Quantum Bit Systems 



In the classical physical world, in classical computation, a bit is represented 
by two states of a physical system (e.g. charge or no charge) |0), In the 
(sub)atomic world, which is governed by the laws of quantum mechanics, we 
have, along with the classical states |0), |1) (such states could be, e.g., spin 
up or spin down of an electron), also their superpositions: 

a |0) + «i ]1) (a , «i G C, |a | 2 + l a i| 2 = 1), 

that is, the linear combination of the classical states |0) and |1). The crucial 
problem however is that if the system is in state «o |0) + a± we cannot 
access, measure «o and cti directly. Instead, a measurement destroys the su- 
perposition, causing the system to return to state |0) with probability |a | 2 
and to state |1) with probability |o;i| 2 . 

The mathematical framework for a quantum bit, that is, a quantum system 
with two classical basis states, is the two dimensional complex Hilbert space 
H\ := C 2 . Let eo,ei be its unit vector basis. Following quantum mechanical 
notation, we write |0) instead of eo and |1) instead of e±. The unit sphere of 
H 1: i.e., the set of all elements of norm 1, is the set of states of the qubit. 

A system of m interacting qubits (m-qubit system) is represented by the tensor 
product 

H m :=H 1 ®H 1 ®...0H 1 , 

S v ' 

m 

which is the 2 m -dimensional complex Hilbert space. It has the canonical basis 

e io <S> e h <g> . . . <g> e im _ x (i , i u . . . , i m _i) G {0, l} m . 
We make further notational conventions: 

e io ® e h (g) . . . (g) e ím _ 1 =: \i ) |* x > . . . \i m -i) =■ \i) , 

where i := (i ii . . . i m _i) 2 := ïfe2 m_1_fe . The vectors \i) (i = 0, . . . , 2 m - 

1) are called classical states, while a general state has the form 

2 m -l /2 m -l \ 

m= E «il<> E N 2 = i ■ 

i=0 \ i=0 ) 

As in the case of a single qubit, measuring an m-qubit system in a super- 
position state \ip) results in one of the classical states \i) with probability 
\ai\ 2 (i = 0, . . . , 2 m — 1). So a state is a linear combination of all possible 
classical states, the coemcients giving the probability that after measurement 
the system moves to this state. 
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One more notational convention: if |-) contains a number inside, or a symbol 
which is used to denote such a number, we mean the canonical basis vector 
corresponding to this number, if |-) contains just a general symbol, like 
we mean any vector of H m (this should be clear from the context). 



1.3 Quantum Computation 



How to use such systems for computing? To make this clear, let us first consider 
an example of a classical computation - the addition of two m-bit numbers, 
which we write as follows: 



\Í ) ■ ■ ■ \im-l) bo) • • • \jm-l) |0> ... |0> 
I 

\Íq) ■ ■ ■ IV-l) |jo) • • • bm-l) |&o) • • • \k m ) 

This computation is realized using circuits of classical gates (and, or, not, 
xor) in the usual way: add the last bits, then the second last plus the carry 
bit etc. Let us emphasize here: Classically, we add two numbers at a time. 

How to operate m-qubit quantum systems? Which operations are allowed? 
Schròdinger's equation implies: all evolutions of a quantum system must be 
represented by unitary transforms of H m . Here is the starting point: 

Quantum computing assumes that we are able to perform a number of ele- 
mentary unitary transforms (quantum gates) on the system. 

What are these operations? Let us consider one Standard set Q m of them. 
First we describe the one-qubit gates - these are gates that manipulate only 
one component of the tensor product H m — H-i® Hi<S> ■ ■ ■ <8> The Walsh- 
Hadamard gate W : H.\ — > H\ is defined by 

(the vàlues on the basis vectors define the unitary transform uniquely). Its 
action on the j-th component of H m is then given by the unitary operator 

= Id® . . . ® Id® ®Id <g> . . . <g> Jd, 

j 

where Id stands for the identity operator on H\. For a real parameter < 6 < 
2n the phase shift Pq : H\ — > H\ is defined as 

P e |0) = |0) P e \l) = e te \l}. 
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We define Pg^ in the respective way. Next we consider two-qubit gates - they 
manipulate any chosen two components of H 1 <g> Hi <g) . . . <g> H 1 . The quantum 
xor gate (also called controlled-not gate) X : H\ ® H\ — * Hi ® H\ is given by 

X|0) |0> = |0> |0) 
X|0)|1) = |0)|1) 
X|1)|0) = |1)|1) 
X|l)|l) = |l) |0) 

That is, if the first bit is zero, nothing happens to the second, and if the 
first is one, the second is negated (controlled not). We can look at this gate 
also as follows: the xor of the two bits replaces the second bit. Denote by 
X^ ,c> : H m — > H m the unitary operator given by applying X to the fc-th and 
í-th. component, that is, 

X m' e) l«o) ■ ■ ■ Vk) ■ ■ ■ \U) ■ ■ ■ \im-l) = No) ...\y) ...\z) ... \Í m -l) , 

where 

\y) \ z ) '■= x \h) \k) ■ 

Now we define 

G m = {wji\P e %XW • 0<j, k ^£<m-l,0<9<2n}. 
The following two results can be found, e.g., in [21]. 

Theorem 1 The set Q m is a universal system of gates - each unitary trans- 
form of H m can be represented as a finite composition of elements of Ç m (up 
to a complex s calar factor). 

The set Q m is still an infinite set. Now consider the following finite subset: 
Sí = {Wg>, P& m , : < j, i < m - 1} . 



Theorem 2 The set forms an approximately universal system of gates - 
each unitary transform of H m can be approximated in the operator norm to 
each precision by a finite composition of elements ofÇ^ (again up to a complex 
factor). 

So once we can implement these gates we can do all unitary transforms (of 
course, the efficiency of such an approximation is still an issue, we will come 
back to that later on). Physicists are working on implementations of these 
gates in various quantum systems such as photons, trapped ions, magnètic 
resonance systems etc. Let me emphasize some crucial points: 
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1. These gates can transform classical states into superpositions. For example, 
the Hadamard gate applied to the first and then to the second qubit: 

|0> |0> — , | ()0> |0) + |0)|1) + |1)|0) + |1)|1)). 

2. They also act on superpositions. For instance: 

2.1 The quantum xor: 

a |0> |0> + a x |0) |1) + a 2 |1) |0) + a 3 |1) |1) 
I 

a |0> |0) + ai |0) |1> + a 2 |1) |1) + « 3 |1> |0) 

2.2 Quantum addition of binary numbers: 

Uij No) • • • \im-l) | jo) • • • \jm-l) |0> ... |0) 

I 

) |jo) • • • lim-l) |^o) • • • |^m) 

T/idí is, in the quantum world, we add all possible binary m-digit numbers 
in parallel (assuming that all ^ 0). So, is a quantum computer an ideal 
parallel computer, with exponentially many processors? Not exactly, it is not 
that easy! We cannot access all components of the superposition. We have to 
measure, thus destroying all results but one! We would get 

\Í0) ■ ■ ■ \Í m -l) |jo) • • • \jm-l) \k ) ■ ■ ■ \Ki) 

with probability | ct^- 1 2 . In this sense a quantum computer is like a black box - 
we can manipulate it, but we cannot "look into it". Anyway, we have reached 
the quantum model of computation (the general way a quantum algorithm 
should look like): 

Quantum model of computation: 

starting state: \i ) G H m , a classical state 
computation: Lq, . . . , U n £ Q m 

\io) -> Ui \i ) -> U 2 Ui \i ) . . . 
:=U n U n - 1 ...U 2 U 1 \i Q ) 
measurement: \ip) = z2i=o a i \V 

— > \i) with probability | q;^ | 2 

repeat 
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The cost of the algorithm is defined to be the number of gates n. We may 
allow for a classical computation after each measurement to determine the 
next starting state, and, after the last measurement, the final result (see [14] 
for a formal approach). Then, speaking about the cost of the algorithm, we 
also have to take into account the needed classical gates. A quantum algorithm 
is said to solve a given problem if it produces a solution with probability > 3/4 
(any fixed constant 7 > 1/2 would do, since the success probability can always 
be increased by suitable repetitions). 

1.4 Computational Power of Quantum Computer s 

What can we do with a quantum computer that we cannot do on a classical 
computer? A simple answer related to Monte Cario algorithms can be given 
as follows: On classical computers, pseudo-random numbers (that is, in fact, 
fully deterministic numbers) are used instead of random numbers. A quantum 
computer would provide us with true randomness. It is simple to let a quan- 
tum computer act as a random number generator: 

Quantum random number generator 

|0)...|0) 

| (Walsh-Hadamard transform) 

i=0 

l (measurement) 
\i) with probability 2~ m 

So, on a quantum computer, we can implement Monte Cario with true (phys- 
ically based) randomness. Let us go one step further: What can be done on a 
quantum computer that a randomized classical computer (that is, a classical 
computer with access to true randomness) cannot do? We already mentioned 
Shor's result [26]: An algorithm which, for a composite integer N, finds a non- 
trivial factor with cost 0((log N) 3 ). No polynomial in logiV classical (deter- 
ministic or randomized) algorithm is known. Recali that there are polynomial 
randomized prime number tests (see Solovay and Strassen [29], and Rabin 
[25]), but they don't give the factors when they find compositeness. 

Although not very likely, it is not excluded that there are polynomial classical 
factoring algorithms (maybe even with cost 0((log N) 3 )). So formally, the 
above result does not strictly prové the superiority of quantum computing. 
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Matters are different with Grover's result [10]. The problem is the following: 
Given / : {0, . . . , N — 1} — > {0, 1} with the property that there is a unique Íq 
with f(io) = 1, find this i . It is not difficult to see that classical (deterministic 
or randomized) algorithms need cost Q(N). Grover [10] produced a quantum 
algorithm that solves the problem with cost 0(y/~N). Here we have to extend 
the above model of computation by what is called a quantum query (quantum 
subroutine, quantum oracle), which is given by the unitary map Qf : H rn — > 
H m defined as 

Qf ■ \i) \y) - \i) \y © f(i)) (o<t< 2 m ~\ y e {o, i}), 

where © is addition modulo 2, and we assumed N = 2 m_1 . Note that if the last 
qubit is set to zero, Qf is like a subroutine which writes the function value 
to the last entry. This is what we essentially assume in the design of most 
numerical and Monte Cario algorithms: the function we want to handle, e.g., 
integrate, is given as a subroutine. (In the quantum setting, addition modulo 
2 guarantees that Qf is a bijection on the classical states and thus, extends to 
a unitary operator.) It can be shown that such a quantum query can always 
be implemented emciently (with only a constant factor increase of resources), 
once we have a classical implementation for /. A quantum computation is 
then still given by the scheme in section 1.3, but with the difference that 
Ui G S m U {Qf}. It is important to mention that the rest of the U should not 
depend on /, that is, the algorithm gets its information about / only through 
query calls. The cost is defined to be the number of query calls, multiplied by 
the number of qubits, plus the number of gates . 



2 Summation and Integration 



Now we consider numerical problems. Let us formulate the general integration 
problem: Given a set D, a measure fi on it, and a fi- integrable function / : 
D — > R, compute (approximate) 

/ f(t)dn(t). 

D 

As its simplest instant, this includes computing the mean (and hence the sum) 
of a finite sequence: 

£> = {0,...,7V-1}, f:D^R, MW) = ^ 
íf(t)dn(t) = ^f{i). 

n *=0 
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Theorem 3 (Brassard, H0yer, Mosca, and Tapp [5]) For all n,N, n < N, 
there is a quantum algorithm which computes the mean j^Y^f=Q l f(i) f or a ^ 
sequences f(i) with f(i) G {0, 1} with cost 0(n) and error 0(1/ n). 

We use the Ò notation to indicate that we suppress possible logarithmic fac- 
tors. The result is easily extended to real valued f(i) G [—1,1]- We have to 
specify how a quantum algorithm can access such vàlues of /. Without loss 
of generality we suppose that N is of the form N = 2 mi . We assume that we 
have a quantum query (subroutine) providing essentially the first m 2 digits of 
f(ï) for some suitable m 2 . More precisely, we can use the unitary mapping Qf 
on H m = H mi <S> H m2 defined as 

Qf ■ 10 \v) -> 10 Iv© [/(0U> (o < % < 2 m \ o< y < 2 m2 ). 

Here © is addition modulo 2 m2 and [f{i)]m 2 is the 7712-bit integer \_2 m2 ~ 1 (f(i) + 
1)J if f(i) < 1 and 2 m2 — 1 if f(i) — 1. (If we have to deal with functions taking 
vàlues in other intervals, say [a,b], we use a similar encoding as a binary 
integer, [f(i)] m2 = [2™*(f(i) - a)/(b - a)\ if f(i) < b and [f(i)] m2 =2 m *-l if 
f(i) = b.) We say that a quantum algorithm computes a number with error e, 
if with probability > 3/4 the result of the algorithm is within e of that number. 
For more details on queries for numerical problems and the error definition 
we refer to [14]. Theorem 3 should be viewed under the aspect of huge iV and 
moderate n. For comparison, let us mention that for classical deterministic 
algorithms the error at cost n for n < N/2 is (even n < cN for any fixed 
< c < 1 suffices). In the classical randomized setting we obtain G(l/y / n) 
(see [16] for details on the latter two statements). 

How to treat integration? Define the Hòlder classes J^ Q for r G N , d G N, 
and < q < 1 as 

r/={fec r ([o,i] d ) : imL<i, 

\d a f(x) - d a f(y)\ < \x - y\e,\a\ = r}. 

Here C r ([0, l] d ) stands for the set of r times continuously differentiable func- 
tions on [0, l] d , || ||oo denotes the supremum norm, a represents a multiindex, 
and d a is the respective partial derivative. 

Theorem 4 (Novak [23]) For alln, there is a quantum algorithm which com- 
putes the integral of functions f G T r ^ e with cost 0(n) and error 
(9( n -(»-+ff)/d-l). 

Because of its relation to Monte Cario methodology it is interesting to describe 
the idea behind the algorithm: We use the technique of separation of the main 
part (in other words, of introducing a control variate) to reduce integration 
to summation in an optimal way Let P n be some interpolation operator with 
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n nodes (tí)™ = q G D = [0, l] d suitable for J^' e in the sense that it gives the 
optimal approximation order 

sup|/(í) - P n f(t)\ = 0(n'^ d ). 
teD 

Represent 

I f(t)dt = I P n f(t)dt + I (f(t) - P n f(t))dt. (1) 

D D D 



Clearly, 



/n— 1 
P n f(t)dt = Y,c l f(n) 

D i=0 



is a quadrature, which can be computed exactly, classically, with 0(n) cost. 
The second part is approximated suitably (see [23] for details) in the form 

1 N ~ l 

j (f(t) - P n f(t))dt « - £ (fiu) - P n f(t t )), 

D i=0 

where N may be much larger than n. Finally we use quantum computation to 
approximate the mean on the right hand side. It can be shown that the order 
of the error is the product of the following contributions: 

n -{r+e)/d ^ n -i 

separation of main part quantum computation 

Note that if the last term of (1) were approximated by Standard Monte Cario, 
we would get the exponent — (r + g)/d — 1/2. 

Compare this result with the classical deterministic setting, where the best 
we can achieve is an error 0(n~^ r+e ^ d ) at cost 0(n) and with the classical 
randomized setting with optimal error 0(n~ ( · r+e ^ d ^ 1 ^ 2 ). Look at the exponents 
when (r + g)/d is small! It follows that the speedup over the deterministic 
setting can be polynomial with arbitrarily large power. Note that in [14], 
Section 6, a class of functions of low smoothness is considered and it is proved 
that for this class the gain of quantum over classical deterministic computation 
is even exponential. A similar situation is observed in [30] for path integration. 

A crucial assumption for the summation algorithm in Theorem 3, and thus also 
for Theorem 4, is the uniform boundedness of the set of sequences (with bounds 
independent of n and N). This raised the problem whether the quantum gain 
was tied to that assumption. In particular, what happens in the case of square 
summability, which is the most important case for Monte Cario, that is, / G 
B2 or / G W2 d (see the definitions below)? Will quantum computation retain 
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its superiority or will Monte Cario catch up? Even more generally, what about 
/ G Bp or / G yVp t d for 1 < p < oo? Let us first look at the discrete problem. 
Denote 

r N-l \ 

The summation problem for this class was settled in [14] (case (i) and (ii)) 
and in [17] (case (iii)). The latter answered a question posed in [16]. 

Theorem 5 Let 1 < p < oo. For all n,N, n < N, there is a quantum 
algorithm which computes the mean ^ Z^Lo* f(i) f or ott sequences f G B^ 
with cost 0(n) and error 

(i) 0{n- 1 ) if 2 < p < oo, 

iii) 0(n- 2+2 / p ) if l<p < 2 andn< VN, 

(iii) O (n~ 2 l p N 2 IP- 1 ) if \<p<2andVN<n<N. 

Again, a comparison to the two classical settings might be illustrative (see [16] 
for details): In the classical deterministic case, for n < N/2, nothing better 
than íl(l) can be obtained, in the classical randomized case the optimal rates 
are 6(n~ 1 / 2 ) if 2 < p < oo, and 6(n~ 1+1 / p ) if 1 < p < 2. Now we consider 
integration in the Sobolev classes d , which are defined by 

W; d = {/a p ([0,l] d ) : ||0°/lk<l, H<r}, 

where r G N, 1 < p < oo, and d a is the weak partial derivative. The following 
result is proved in [15]. It answers another question from [16]. The approach 
consists of a new discretization technique, by which one can derive (optimal) 
integration algorithms for Wp d from (optimal) summation algorithms for B^ . 

Theorem 6 Let 1 < p < oo, r, d G N 7 r/d > 1/p (Sobolev embedding condi- 
tion). For all n, there is a quantum algorithm which computes the integral of 
functions f G Wp d with cost 0(n) and error 0(n~ r ^ d ~ 1 ). 

In the classical deterministic setting the optimal rate is 0(n~ r l d ), while in 
the classical randomized setting we have 0(n~ r t d ~ 1 / 2 ) if 2 < p < oo, and 
Q^ n -r/d-i+i/p^ jf i < p < 2 (see [16], [13], and the references therein). The 
quantum rate for 1 < p < 2 comes as a surprise. After previous results one was 
tempted to conjecture that the quantum setting could reduce the exponent 
of the classical randomized setting by at most 1/2. Now we see (p — 1) that 
there can even be a reduction by 1. 

Are these results about quantum algorithms optimal? In other words, is it 
possible to improve the rates by other, better quantum algorithms? To verify 
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optimality, we have to establish lower bounds vàlid for all possible quantum 
algorithms. It turns out that all the results about summation and integration 
presented here are optimal (up to logarithmic factors, at least). The following 
was the first matching lower bound for summation. It shows that the rate in 
Theorem 3 is optimal. 

Theorem 7 (Nayak and Wu [20]) There are constants C\,C2 > such that 
for all n, N with n < C\N the following holds: Each quantum algorithm which 
computes the mean X^Lq 1 /(*) f or a ^ sequences f(i) with \ f(i)\ < 1 using 
at most n quantum queries has error not smaller than Cijn. 

Using their technique and methods of information-based complexity theory, 
the following can be shown (see [23], [14], [17], [15]). 

Theorem 8 The rates established in Theorem 4 (integration in J~2 e )> Theo- 
rem 5 (mean of sequences in B^ ), and Theorem 6 (integration in Wp d ) are 
optimal (up to logarithmic factors). 



In the following table we summarize the results. The respective entries give 
the optimal rates at cost n, constants and logarithmic factors are suppressed. 
The constant < c < 1 in the first column does not depend on n and N. 





deterministic 


random 


quantum 


B^,2<p<oo,n<cN 


1 


n -l/2 




B*, 1 < p < 2, n < y/N 


1 


n -l+l/p 


n -2+2/p 


B£, 1 < p < 2, y/N < n < cN 


1 




n -2/p N 2/p-l 


B?,n<y/N 


1 


1 


1 


B?,VN <n<cN 


1 


1 


n~ 2 N 


S d 


n -(r+g)/d 


n -(r+g)/d-l/2 


n -(r+g)/d-l 


W; d , 2 < p < oo 


n -r/d 


n -r/d-l/2 


n -r/d-l 


w; id , i< v < 2 


n -r/d 


n —r/d-l+l/p 


n —r/d—l 




n -r/d 


n -r/d 


n -r/d-l 
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